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C ' ABSTRACT 

c3 : 
. 

■rj- \ Near-infrared and optical imaging of HD32297 indicate that it has an edge- 

on debris disk, similar to (3 Pic. I present high resolution optical spectra of the 
Na I doublet toward HD32297 and stars in close angular proximity. A circum- 
\ stellar absorption component is clearly observed toward HD32297 at the stel- 

lar radial velocity, which is not observed toward any of its neighbors, including 
the nearest only 0'9 away. An interstellar component is detected in all stars 
>90 pc, including HD32297, likely due to the interstellar material at the bound- 
ary of the Local Bubble. Radial velocity measurements of the nearest neighbors, 
Q-V BD+07 777s and BD+07 778, indicate that they are unlikely to be physically 

associated with HD32297. The measured circumstellar column density around 
\ HD32997, logA^Nai ~ 11.4, is the strongest Na I absorption measured toward 

any nearby main sequence debris disk, even the prototypical edge-on debris disk, 
(3 Pic. Assuming that the morphology and abundances of the gas component 
around HD32297 are similar to (3 Pic, I estimate an upper limit to the gas mass 
c3 ' in the circumstellar disk surrounding HD32297 of ~0.3 M ffi . 

Subject headings: circumstellar matter — ISM: structure — line: profiles — 
planetary systems: protoplanetary disks — stars: early-type — stars: individual 
(HD32297) 



1. Introduction 

Debris disk systems provide a look at an intermediate stage of stellar system evolution. 
They represent the transition between the early formation of stars and planets in a primordial 
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protoplanetary disk as seen toward pre-main sequence stars, and the mature stage of an 
evolved system, like our solar system, which is clear of all primordial material and retains 
only a hint of secondary products (e.g., zodiacal dust), the final remnants of the stellar and 
planetary formation process. Although a debris disk has lost most of its primordial material, 
the observed infrared luminosity of circumstellar dust, caused by collisions of planetismals 
and other small bodies, is typically several order s of magnitude larger than estimated for the 
Kuiper and asteroid belts in our solar system (IBackman fc Parescd I1993T ). Ever since the 
detection of dusty ci rcumstellar material around main sequence stars via infrared excesses 
( Aumann et al.lll984T) . researchers have been looking for circumstellar gas phase absorption 
( IHobbs et al.lll985l ). Of the initial major infrared excess main sequence stars, only Pic 
showed gas phase absorption in optic al absorption lines (e.g ., Ca II and Na I), due to its disk 
morphology and edge-on orientation (ISmith fc Terrildll984j ) . Such on orientation provides a 
unique opportunity to simultaneously measure both the dust and gas components of a debris 
disk, at an interesting transition near the end of stellar and planetary formation. 



Only a few oth er edge-on debris di sks have been found since, including C ar (ILagrange-Henri et al 



1990 bh. HD85905 flWelsh 



et al 



19981 . HR10 flLagrange-Henri et al.lll990ah. and AIJ Mic 



fcalas. Liu, fc Matthewslbooi iRoberge et alil2005h . iRedfield. Kessler-Silacci. fc Cieza] ^0071 ) 
observed Car, HD85905, HR10 with the Spitzer Space Telescope and did not find strong 
infrared excesses toward any of them, although an optical monitoring campaign showed 
clear signs of gas variability, as noted by researchers earlier. However, the magnitude of 
circumstellar absorption in these systems is lower than observed toward Pic. 



Long Ca II monitoring campaigns of Pic (e.g., iPetterson fc Tobinlll999l ). find signifi- 
cant short-term absorption variability. This variability can be explained by gas clouds very 
close to the star, which are caused by evaporatin g, star-grazi ng, km-sized objects, simply re- 
ferred to as, Falling Evaporating Bodies (FEB's; lBeustlll994j ). A strong "sta ble" component, 
at re st in the stellar reference frame, is also detected toward Pic (e.g., ICrawford et al. 
19941 ) . The distribution of gas in this component, contrary to the variabl e component lo 
cated very close to the star, is dispersed throughout the extended dust disk (IBrandeker et al. 
2004h . 



A "stable" absorption component in a gas phase resonance line can be caused by ei- 
ther intervening circumstellar or interstellar gas. Measuring the interstellar medium (ISM) 
along the line of sight and in the locality surrounding a circu mstellar disk candidate, is crit 
icalto characterizing any "contaminating" ISM absorption (I Crawford! 1200 ll ; IRedfield et al. 
20071 ). In particular, the Sun resides in a large scale ISM structure known as the Local Bub- 
ble, whose boundary a t ~100pc is defined by a significant quantity of interstellar material 
(jLallement et al.ll2003l ). If a "stable" absorption component is observed at the stellar radial 
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velocity, and similar absorption is not detected toward any proximate stars, it is likely that 
the absorption component is caused by circumstellar material. 

Using near-infrared scattered ligh t observations taken with the Hubble Space Telescope, 
Schneider. Silver stone, fc Hined (120051 ) discovered that the debris disk surrounding HD32297 
has an edge-on orientation. Disk emission extends out to ~400AU in their observations, 
while radii <33. 6AU are occ ulted by the coronagraphic obstacle. Optical scattered light 
observations by iKalad (120051 ) confirmed this orientation and extended the range of disk 
emission to ~1680 AU. The edge-on orientation of HD32297 makes it an ideal target for gas 
phase absorption measurements. 



2. Observations 



Observations of the Na I D doublet (5895.9242 and 5889.9510 A) toward HD32297 were 
made over several epochs. The Na I doublet is among the strongest trans itions in the optical 
wavelength band, appropriat e for observing interstellar faedfieldl 120061 ) and circumstellar 
( ILagrange-Henri et al.l Il990bl ) absorption toward nearby stars. In addition, several stars 
in close angular proximity to HD32297 were observed, in order to reconstruct the ISM 
absorption profile along the line of sight. Stellar parameters of the observed targets are 
given in Table dj an d the observational parameters are listed in Table [2J 

High resolution optical spectra were obtained using the Coude Spectrometer on the 
2.7m Harlan J. Smith Telescope at McDonald Observatory. The spectr a were obtained at 
a resolution of R = A/AA ~ 240,000, using the 2dcoude Spectrograph (jTull et al.lll9951 ) in 
the CS2 1 configura tion. The data were reduced using Image Reduction and Analysis Facility 
(IRAF; lTodylll993l ) and Interactive Data Language (IDL) routines to subtract the bias, flat 
field the images, remove scattered light and cosmic ray contamination, extract the echelle 
orders, calibrate the wavelength solution, and convert to heliocentric velocities. Wavelength 
calibration images were taken using a Th-Ar hollow cathode before and after each target. 

Numerous weak water vapor lines are commonly present in spectra around the Na I 
doublet, and must be modeled and removed, in order to measure an accurate interstellar (or 
circumstellar) Na I abs orption profile. I use a forward modeling technique demonstrated by 
Lallement et al.l (119931 ) to remove telluric line contamination in the vicinity of the Na I D 
lines, with a terrestrial atmosphere model (AT - Atmospheric Transmission program, from 
Airhead Software, Boulder, CO) developed by Erich Grossman. With two Na I absorption 
lines, it is straightforward to identify contaminating telluric absorption. 



All absorption lines were fit using standard methods (e.g., §2.2 in iRedfield fc Linsky 
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2004J). Gaussian a bsorptio n comp onents are fit to both Na I D lines simultaneously using 



atomic data from iMortonl (120031 ). and then convolved with the instrumental line spread 
function. Fitting the lines simultaneously reduces the influence of systematic errors, such as 
continuum placement and contamination by weak telluric features. The free parameters are 
the central velocity (v), the line width or Doppler parameter (b), and the column density 
(N) of Na I ions along the line of sight. The fits are shown in Figured] and fit parameters 
with la statistical errors are listed in Table [2j 

In addition, the spectra were used to estimate the stellar radial velocity (i>r) and pro- 
jected stellar rotation {vsini) for HD32297, BD+07 777s, and BD+07 778 (see Tabled]), 
quantities not listed in SIMBAD for these targets. The radial velocities of all 3 objects differ 
significantly, and therefore it is unlikely that they are physically associated. Note that the 
radial velocity of HD32297 (t>R ~ +20 km s -1 ) is measured from broad Na I and Ha stellar 
absorption lines, and therefore is not tightly constrained. 



3. Identification of Circumstellar Absorption Toward HD32297 



The left column of Figure [T] shows that Na I absorption is clearly detected toward 
HD32297 in 5 observations over 5 months. Two components are easily distinguished, a strong 
component at ~24.5 km s _1 and a weaker component at ~20.5 km s" 1 . The Na I spectral 
region for 5 stars in close angular proximity to HD32297 is also shown in Figure [H Only a 
single ISM component, at ~24.2 km s _1 , is detected in the 3 distant neighbors, indicating 
that large scale interstellar material is located at a distance between 59.4-112 pc. All targets 
located beyond this material, including HD32297, should have a similar ISM absorption 
feature. This strong ISM absorption is probably associated with the bo undary material 
of th e Local Bubble, which is estimated to be ~90 pc in this direction (jLallement et al. 



20031 ). If located at this distance, the physical separation of the interstellar material observed 
toward HD32297 and the material toward BD+07 777s (AO = 0'.9) is 0.025 pc, BD+07 778 
(2'A) is 0.064 pc, and 18 Ori (5?1) is 8.1 pc. Toward HD32297's two closest neighbors, the 
ISM absorption is almost identical in projected velocity and column density to the strong 
absorption seen toward HD32297, while toward 18 Ori, the absorption differs slightly in both 
v and N, indicating that any small scale morphological variations in the Local Bubble shell 
are on scale s >0.1pc but <8pc. Sma l l scale variations in the Local Bubble shell have been 
detected by lRedfield. Scalo. fc Smith] (120061 ) on scales ~0.5pc. 



It is unlikely that the unique 20.5 km s 1 feature observed toward HD32297 is caused by 
a small scale i nterstellar structure. Although sma l l ISM structures (0.01- 2.0 pc) have been 



observed (e.g., iFerlet. Dennefeld. Maurice! Il985l ; iMeyer & Blades! Il996l ). it is more likely 
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that the unique feature is due to absorption in the circumstellar environment surrounding 
HD32297 because (1) HD32297 is known to be an edge-on debris disk, (2) no similar absorp- 
tion is detected in the very close neighboring sightlines (0.03-0.08 pc), and (3) the absorption 
matches the stellar radial velocity. 



4. Temporal Variability of Circumstellar Component 



Temporal variability is also a hallmark of circum stellar material (e.g-. lFerlet. Vidal-Madjar. fc Hobbs 



19871 . IPetterson fc Tobinlll999l . iRedfield et al.ll20071 ). To search for variability, Figure [2] shows 
difference spectra of all observations. Some indication of temporal variability on time scales 
of months is detected. For example, between the 2005 Sep and 2006 Feb observations, the 
~20.5 km s _1 feature became stronger and the separation between the circumstellar and in- 
terstellar features became less distinct, despite the fact that the 2006 Feb observations were 
made at a slightly higher resolving power. The redshifted variability seen between 2005 Sep 
and 2006 Feb is ~5<j above the standard deviation. The same pattern is seen in both Na I 
lines, indicating that the telluric contamination is not causing the variation. Slight changes 
in the resolving power of our instrument could mimic this variable behavior, differentially 
moving light from the cores of the line to the wings (or vice versa). However, resolution vari- 
ability should cause (1) symmetric features in the wings of the line whereas we see a feature 
only to the blue of the ISM feature and not to the red, and (2) should have a stronger effect 
on stronger absorption features, whereas the feature is roughly identical in both lines, which 
could be caused if the a bsorbing material covers only a fraction of the stellar disk, as has 
been seen toward j3 Pic ( 1 Vidal-Madjar et al.lll994l ). 



This data alone provides only a subtle indication of temporal variation in Na I, partially 
because any significant absorption toward the red, is masked by the strong ISM feature. 
Redshifte d circumstellar absorptio n dominates the Ca II gas absorption variability toward j3 



Pic (e.g., IPetterson fc Tobinlll999l ). while no temporal variability has ever bee n detected in 



Na I toward (3 Pic, only the "stable" absorption component is seen in this ion (IWelsh et al. 



19971 ). Circumstellar variability in Na I has been detected in other edg e-on debris disks, 
e.g., /3 Car, HD85905, and HR10 JWelsh et al.lll998l : IRedfield et all 120071 1 . Any redshifted 
absorption occuring in this object could cause fluctuations in the measured column density 
of the "constant" ISM feature. Little evidence for variability is found toward the blue. 
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Gas Disk Mass 



These observations indicate that HD32297 has the strongest Na I circumstellar disk sig- 
nature detected around a nearby main sequence debris disk star. Even compared to j3 Pic, the 
protot ypical edge-on debris disk with Na I absorption column densities of logiV NaI ~ 10.69- 
10.73 flHobbs et al1ll985l : IvVelsh et al.lll997h . the gas disk around HD32297, with log A^ NaI ~ 
11.4, has 5x the Na I column density. A crude estimate of the gas mass surrounding HD32297 
can be made if it is assumed to have the same morphology and abundances as the stable gas 
around (3 Pic. Although the observations of HD32297 indicate some red-shifted temporal 
variability, much of the gas is stable over all observations. Using Pic as a proxy, the y ariable 
gas is likely located very close to the star (lLagrange. Backman. fc Artvmowiczll2000h. while 



the s table gas at rest in the stellar frame, likely traces the bulk dust disk (IBrandeker et al. 



20041 ) . For this calculation, I assume all the gas is in the stable component, and there- 
fore this gas mass estimate should be considered an upper limit. The morphology of the 
disk is assumed to follow a broken power la w density profile, as fit to the Na I emission 
profile of the f3 Pic disk (see Equation 1 of IBrandeker et all I2004T) . and assumed to ex- 
tend out to the edge of the debris disk at ~1680AU kalasll2005h. The abundances in the 
HD32297 disk are assumed to be similar to (3 Pic (IRoberge et al.l 120061 ), where the ratio 



N(H 1)/N(N& I) < 8.8 x 10 8 , is based on P ic Na I measurements by IBrandeker et al. 



(120041 ) and H I limits by iFreudling et al.l (119951 ). Given these assumption, I calculate a gas 
mass, distributed through the bulk debris disk surrounding HD32297 at M gas ~ 0.3M e . 

Future observations are planned to continue monitoring the temporal variability of the 
circumstellar gas toward HD32297 to determine the ratio of stable to variable gas, and 
measure the Ca II gas disk absorption, in order to independently measure the Ca II to Na I 
ratio. A more definitive detection of temporal variability may require monitoring excited 
lines which will show circumstellar absorption, but not the strong interstellar feature. 



6. Conclusions 

I present the first high resolution optical spectra of the Na I doublet toward the debris 
disk HD32297 and stars in close angular proximity. A summary of results include: 

(1) Two absorption components are detected toward HD32297, while only one is detected 
toward its proximate neighbors located at a comparable distance. The extra absorption com- 
ponent in the spectrum of HD32297, which is also at rest in the stellar reference frame, is 
therefore likely caused by circumstellar material. 

(2) The ISM absorption is similar among HD32297 and its two closest neighbors, and is 
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likely due to absorption from the shell that defines the boundary of the Local Bubble. Some 
variation in Local Bubble absorption is detected toward 18 Ori. 

(3) Radial velocities of HD32297, BD+07 777s, and BD+07 778 are measured and differ 
significantly, indicating that they are likely not physically associated. 

(4) Some indication of temporal variability is detected over several epochs of observa- 
tions. Instrumental resolution variations and masking by the strong ISM absorption, make 
a definitive detection of circumstellar Na I variability difficult. 

(5) The measured circumstellar feature toward HD32297 (log iV NaI ~ 11.4) is the strongest 
such absorption measured toward any nearby main sequence debris disk, ~5 times greater 
than the column density of the prototypical edge-on debris disk, (3 Pic. 

(6) If the morphology and abundances of the stable gas component around HD32297 are 
assumed to be similar to f3 Pic, I estimate an upper limit to the gas mass in the circumstellar 
disk surrounding HD32297 of ~0.3 M e . 

Support for this work was provided by NASA through Hubble Fellowship grant HST- 
HF-01 190.01 awarded by the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 
5-26555. I would like to thank D. Doss, G. Harper, and A. Brown, for their assistance with 
these observations. The insightful comments by the anonymous referee were very helpful. 
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HD32297 Neighbors 




Na I A5889 Velocity (km s" 1 ) Na I X5889 Velocity (km s" 1 ) 

Fig. 1. — Na I absorption lines toward HD32297 (left) and toward stars proximate to 
HD32297 (right). All spectra are flux normalized such that stellar features are removed, 
except for BD+07 778, where the stellar feature is retained. All objects beyond 60 pc, in- 
cluding HD32297, show a significant interstellar component at ~24.5 km s _1 , whereas only 
HD32297 has an additional absorption feature ~20.5 km s _1 , presumably circumstellar. 
Temporal variability can be seen in the circumstellar component and redward in HD32297 
over several observational epochs. The absorption in both lines of the Na I doublet is fit 
simultaneously. The best fit, convolved with the line spread function, is overplotted (red). 
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Fig. 2. — Difference spectra of all HD32297 observations for both Na I lines, relative to the 
final spectrum taken 2006 Feb 15, which is plotted as the dashed line. The gray shaded region 
indicates the ratio error band. The ratio of the telluric water vapor spectra of 2005 Sep 15 
and 2006 Feb 15 is shown offset toward the bottom of each plot. The most notable temporal 
variation occurs between the 2005 Sep 15 observation (red), where the flux variation, centered 
at ~22.5 km s _1 , is ~5cr. This redshifted variability is common in gas absorption toward 
edge-on debris disks, although in this particular instance the variability is partially masked 
by the strong ISM component. 
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Table 1. Stellar Parameters for HD32297 and Nearby Stars 1 



HD 


Other 


Spectral 


my 


vji v sin i I 


b 


Distance 13 


A6» 


^^*pos C 


# 


Name 


Type 


(mag) 


(km s _1 ) (km s -1 ) (deg) 


(deg) 


(pc) 


(deg) 


(PC) 



32297 BD+07 777 AO 8.13 ~+20 d ~80 d 192.83 -20.17 112±i| 0.0000 0.0 





BD+07 777s 


G0 C 


10.2 


-55 d 


~2 d 


192.85 


-20.17 




0.0156 


0.030 


32304 


BD+07 778 


G5 


6.87 


-1.4 d 


~3.5 d 


192.88 


-20.17 


13413 


0.0406 


0.079 


30739 


7T 2 Ori 


AlVn 


4.35 


+24 


212 


189.82 


-21.83 


59.41" 


3.2653 


3.4 


31295 


7T 1 Ori 


A0V 


4.66 


+13 


120 


189.35 


-20.25 


37.0+J-^ 


3.2748 


2.1 


34203 


18 Ori 


A0V 


5.52 


-8.2 


70 


191.29 


-15.25 


112.9l^ 8 4 


5.1306 


10.1 



a All values from SIMBAD unless otherwise noted. 
b Distances calculated from Hipparcos parallaxes. 

c Physical separation in the plane of the sky from HD32297, at the distance of the closest partner. 
d Measured from spectra presented in this paper. 
c Spectral type based on B — V = 0.57. 



Table 2. Observational and Absorption Fit Parameters for HD32297 and Nearby Stars 



HD 


Other 


Date 


fatm a 


(FWHM ThAl ) h 


S/N c 


V 




b 




log N 


# 


Name 




(km s- 1 ) 


(km s- 1 ) 




(km s~ 




(km s 




(cm- 2 ) 


32297 


BD+07 777 


2005 Sep 15 


+28.2 


1.565 


31 


20.459 ± 
24.546 ± 


0.030 
0.012 


0.38 ± 
0.742 ± 


0.28 
0.063 


10.97i°; 2 ? 

1 O QQO~T"0.083 

12.dS8_ 103 


32297 


BD+07 777 


2006 Jan 26 


-22.8 


1.253 


17 


20.50 ± 


0.27 


1.23 ± 


0.55 


11.384 ± 0.077 














24.30 ± 


0.12 


0.84 ± 


0.21 


12 16+ 012 

11 44O+0.053 


32297 


BD+07 777 


2006 Jan 28 


-23.4 


1.266 


22 


20.48 ± 
24.34 ± 


0.23 
0.17 


1.24 ± 
0.82 ± 


0.42 
0.31 


32297 


BD+07 777 


2006 Jan 29 


-23.7 


1.288 


20 


20.42 ± 
24.424 ± 


0.20 
0.086 


1.71 ± 
0.71 ± 


0.28 
0.23 


11 454+0-049 

L1.101_ Q 050 

1 2 D2 +0 ' 18 


32297 


BD+07 777 


2006 Feb 15 


-27.5 


1.199 


39 


20.51 ± 
24.41 ± 


0.26 
0.13 


1.14 ± 
0.94 ± 


0.41 
0.30 


11 297+ - 097 
-0.120 

11 97+0- 15 





BD+07 777s 


2006 Feb 16 


-27.7 


1.238 


6 


25.55 ± 0.28 


1.581 ± 0.078 


i2.03^;; : j 2 


32304 


BD+07 778 


2006 Feb 16 


-27.7 


1.244 


52 


24.583 ± 0.083 


1.65 ± 0.15 


i2-i8±g:ii 


30739 


7T 2 Ori 


2006 Feb 16 


-28.3 


1.213 


133 






< 9.8 


31295 


tt 1 Ori 


2004 Oct 18 


+21.1 


1.761 


100 






< 10.1 


34203 


18 Ori 


2006 Feb 15 


-27.3 


1.189 


117 


22.464 ± 0.022 


1.150 ± 0.036 


11.4U.5_q 044 



a Projected velocity of the Earth's atmosphere. 

b Resolution based on neighboring ThAr comparison lamp spectra assuming that the ThAr lines are fully resolved. 
c Signal-to- noise at core of Nal absorption line. 



